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ABSTRACT 

Using the most recent proper- mot ion determination of the old, Solar- 
metallicity, Galactic open cluster M 67, in orbital computations in a non- 
axisymmetric model of the Milky Way, including a bar and 3D spiral arms, we 
explore the possibility that the Sun once belonged to this cluster. We have per- 
formed Monte Carlo numerical simulations to generate the present-day orbital 
conditions of the Sun and M 67, and all the parameters in the Galactic model. 
We compute 3.5 x 10 5 pairs of orbits Sun-M 67 looking for close encounters in 
the past with a minimum distance approach within the tidal radius of M 67. In 
these encounters we find that the relative velocity between the Sun and M 67 
is larger than 20 km/s. If the Sun had been ejected from M 67 with this high 
velocity by means of a three-body encounter, this interaction would destroy an 
initial circumstellar disk around the Sun, or disperse its already formed plan- 
ets. We also find a very low probability, much less than 10~ 7 , that the Sun was 
ejected from M 67 by an encounter of this cluster with a giant molecular cloud. 
This study also excludes the possibility that the Sun and M 67 were born in 
the same molecular cloud. Our dynamical results convincingly demonstrate that 
M67 could not have been the birth cluster of our Solar System. 

Subject headings: (Galaxy): open clusters and associations: individual: NGC 2682 
(Messier 67) — Galaxy: structure 
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Introduction 



In addition to its intrinsic interest for being the birth place of the only life form we know, 
the origin of the Solar System provides important constraints to the current paradigms of 
star and planet formation. Identifying some of the history of the Earth's climate (ice ages, 
extinctions, etc.) to our place in the Galaxy, knowing some of the Sun's orbital dynamic 
history would shed some light to important matters on different disciplines, not only Astron- 
omy. In this direction, an interesting question is how the stellar birth environment influences 
the formation and evolution of a planetary system. 

Although there is a fraction of stars that was likely formed in isolated environments, 
studies f rom the last decade s have revealed that most stars were born within groups and 



studies i rom tne last decade s nave revealed tnat most stars were born witnm group ; 
clusters (Eada fc Ladal Il995l: ICarpenterl I2OO0I : kdams fc Mversl l200ll ; lLada fc Ladal 



Porras et al 



2003; 



Allen et al. 



2003 



20071 ). and the Sun does not seem to be the exception. 



There are several indications that the Solar System was formed in a gravitationally bound 
cluster; these inclu de the observed ex citation of the Kuiper belt and the extreme orbital 
elements of Sedna (IBrown et al.l I2004J ) and some other bodies in the Kuiper belt. More- 
over, the presence of short-lived radioisotop es in meteorites indicates t hat the Solar System 
formed close to at least one massive star (IWadhwa fc Russell! |2000| ; iGoswami fc Vanhala 
20001 ; iMever fc ClavtorJ I2OO0I : Irlester et all 120041 lAdamsl boioh . Several papers have been 
devoted to investigate the size of this cluster based on the observational constraints posed 
by our planetary disk, such as the Kuiper belt, which represents a very fragile and restrictive 
entity. 

In a relatively dense star cluster, with a stellar density « 1000 pc~ 3 , the Solar System 
shoul d be able to survive for as long as 250 Myr before its disruption becomes likely ((Adams 
20101 ). which can also be inferred from the fact that the incli nation angles of the orbits of 
Neptune and Uranus are not due to perturbations by flybys (IGaidoa Il995l ) . Thus, on one 
hand, considering the long timescales involved for the Sun to have a mild encounter, it seems 
necessary th at it was form ed in a long-lived cluster -which occurs for about 10% of stellar 
population- ( lAdamsl I2010I ). Only relat ively large bound c l usters (of more than a thousand 
stars), are expected to live that long (IKroupa et al.l 1200 It lLamers et al.ll2005|). Therefore. 
the Solar System could have been born in a very large cluster and survive~ lscally fc Clarke 
200ll ; iMann &; Williamsll20091) . provided, for example, that it spent enough time in the outer 
regions (lAdamsl |2010| ). On the other hand, an encounter that took place more than 10 
Myr after the Oort cloud started forming, which probably would haye occurred in a dense 
cluster, would have compromised its formation (ILevison et al.l 120041 ; lAdamsl |201Q| ). In the 
past several years, observations have started to place constraints on the plausible distance 
ranges to clusters in the Solar neighborhood that could have give birth to the Sun (those 
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within 2 kpc from the Sun). If our parent cluster was relatively small, as others in the 
Sol ar neighborhood, it probably contained a few thousand stars, and had a size of a few 



pc (lAdams fc Laugh 



..200 1; Ada msll2010l ). Thus, inferred from the need of a likely nearby 
Adamsl ( 2010 ) proposes that the Sun was probably not too far from the 



supernova explosion, 

cluster center ( ~0.2 pc ) , whe re the density is higher and the more massive stars are placed. 



In his Table 2, lAdamsl (120101 ) summarizes some Solar System properties and constraints on 
the Solar birth cluster. 

The history of the Sun (and its siblings) in the Galaxy, from its birth to the present has 
long been a subject of general interest. In this direction, identifying some of the siblings of 
the Sun would provide a direct observational constraint on its birth cluster, with constraints 
to the number of stars in the cluster, and even to the IMF i f siblings were found over a ra nge 
of stellar masses ( IPortegies Zwartll2009l ; iBrown et al.ll2010l ; iBland-Hawthorn et al.ll2010f ). 



In his work, IPortegies Zwartl ( 120091 ) considered the constraints of the Sun's birth cluster 
and simulated the orbital evolution of stars in a 1-pc virial-radius dissolved cluster, with 
total mass of ~ 10 3 M , along a Sun-like orbit in an axisymmetric potential model for the 
Galaxy. He concluded that, depending on how quickly the cluster became unbound, between 
10% and 40% of th e Sun's siblings should be located within 1 kpc of its present location. 
Brown et al.l ( 120101 ) simulated the orbits in the Sun's birth cluster, starting from an assumed 
birth place for the Sun, obtained from tracing back the Sun's orbit over 4600 Myr in a 
simplified axisymmetric Galactic potential. They generated a birth cluster and integrated 
forward in time in order to find the present day phase space distribution of the siblings that 
remain close to the Sun. Using this phase space distribution, they made a first attempt to 
identify candidate siblings of the Sun by searching in the Hipparcos Catalog ue, and in the 
Geneva-Copenhagen survey of the Solar neighborhood (IHolmberg et al.ll2009l ). They did not 
find convincing Solar siblings within 100 pc from the Sun, although they only examined a 
small fraction of the nearby stars. 



Making a two dimensional, 2D, numerical experiment like that done by IPortegies Zwart 



(2009), but now using a Galactic potential with stationary and transient spiral arms, iMishurov fc Acharova 

( 1201 ll ) show that due to the effect of the spiral arms there is an appreciable drift of the stars 

from the original position of the parent cluster, and only if this cluster had ~ 10 4 stars there 

is a good chance of findi ng Solar siblin g s wit hin 100 pc from the Sun. Also including a 

2D spiral arms potential, iBobylev et al.l (120111 ) have found two possible Solar siblings in a 

sample of 162 candidate stars from the Hipparcos Catalogue, showing that the spiral arms 

favor the close encounters of the Sun with these two stars. 



Now, an interesting question is: if we come from a stellar cluster, is part of this cluster 
still there? Starting from this idea and considering the very similar metallicities, ages, and 
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distance from the Galactic center of the Sun and M 67, it becomes tempting to place the 
Sun origin within this open cluster. 

One of the peculiar and perhaps most intriguing aspects of M 67 is its chemical com- 
position similar to that of the Sun. All recent high quality works, based on observations of 
both, evolved and Solar type stars, indicate an impressive similarity between the M 67 chem 



ical composition and t he Solar one (jTautvaisiene et al.ll2000t iRandich et al.ll2006l ; iPace et al. 



20081 ; iFriel et al.ll2010l ). The similarity in composition and age is so close that several authors 



have i ndicated in M 67 st a rs the best-ever Sol ar analogues so far discovered (IPasquini et al. 



20081 ; lOnehag et al.l l2011t ICastro et al.l 120111 ). closer than any other field star in the So- 
lar neighborhood. Such a similarity is particularly interesting in the context of this work, 
because it is very tempting to apply to the Sun the concept of chemical tagging, and to 
associate therefore our star to M 67. 

Chemical tagging postulates, on the basis of the similarity of chemical composition 
observed in open clusters and stella r groups, that chemical composition can be used to de- 
termine the common origin o f stars ([Freeman fc Bland-Hawthornl|2002t iDe Silva et al.ll2007l ; 
Bland-Hawthorn et al.l 120101 ). This very interesting idea is supported by the evidence that 
most stars in open clusters share the same chemical composition, within th e observational 
uncertainties ( IPasquini et al.ll2004j ; IRandich et al.ll2006l ; IDe Silva et al.l 120071 ). Most impor- 
tantly, this is not limited to Fe or metals, but is also valid for neutron capture elements. 
Similarity in S and R processes are very important because these elements are sensitive to 
the local environment, therefore to the specific location of the star birth. The cluster chemi- 
cal homogeneity does not yet demonstrate that if two stars share the same abundance, they 
share the same birthplace. However , by compiling chemi cal abundances for many nearby 
open clusters and comparing them, IDe Silva et al.l (120091 ) show that large differences are 
present in the elemental abundances of these clusters, strongly indicating that the variabil- 
ity of elements such as Mg, Na, Zr and Ba, appears quite high even for clusters within ~1 
kpc or less from the Sun. Even if the spread of the data can be influenced by the compilation 
of different authors, and by possible systematic trends between computations of the abun- 
dances in dwarfs and giants, the difference seems large enough to guarantee that clusters 
with similar [Fe/H] abundances may well have different chemical composition, especially as 
far as neutron capture elements are concerned. This would suggest that the probability for 
two stars to share the abundance pattern should be low if they were not born from the same 
cloud. 

There is no doubt that, according to this criterion, M 67 is by far the best known 
birth-place candidate for the Sun. It has been known for a long time that the Sun is 
slightly more metal rich with respect to the majority of the stars in the Solar neighborhood 
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( IWyse Gilmord Il995l see however Casagrande et al. 2011, for a different result) and 
the hypothesis that the Sun has moved from its birth radius has been studied in the past 
dWielen et al1ll996l ). 



Respecting the age similarit y between the Sun and M 67, the age of the Sun has been de 



term i ned in several studies (e.g . . iGuenther 



1997 



Dziembows 



1989 



a et al. 1999; Bonanno et al. 



Bahcall et al.l 



2002 



1995 



Baker et al. 



Guenther fc Demarque 



2005 



Moynier et al. 



2007 



Houdek fc Goughll201ll ). Its mean value is around 4.57 Gyr, with a small uncertainty (« 
1 Myr) for the formation of the Solar System, but f or the Sun itself thi s unce rtainty be- 
ing « 0.1 Gyr, or even larger, 0.2 Gyr, according to iDziembowski et al.l (1 199 91). The age 
of M 67 has also be en estimated, with several results: 



Gvr (ISaraiedini et al. 



( VandenBerg fc Stetson 



( iDinescu et al. 



1995 



201 lh : 3.871^ Gyr (ICastro et al.ll201lh : 3.9 ± 0.1 Gyr (Bellini et al. 



e^ 



4.2 ± 



2009); within 3.5 and 4.8 Gyr (lYadav et al 



20041 ); wit hin 3.7 and 4.5 Gvr (Richer et al. 



40 +i.o 



2010aD ; within 3.5 and 
l2008h : 4.0 ± 0.4 Gyr 
3; 4.0 ± 0.5 Gyr 

„ 05 Gyr (jDemarque et al.lll992l ). The mean age of M 67 is ~ 4.0 



.6 Gy r ( IQnehag et al. 



Gyr with a mean uncertainty f^O.5 Gyr. Thus, the Sun's age is approximately within the 
2a error bar of M 67's age, and with their similarity in metallicity, there is some probability, 
to be quantified, that the Sun may have born in M 67, or close to this cluster. 

In this work we investigate the reliability of this hypothesis by computing backward in 
time the orbits of the Sun and M 67 in a Milky- Way-like Galactic potential, including the 
effects of the spiral arms and a Galactic bar, looking for close encounters. Our study considers 
a three dimensional potential for the spiral arms, with the orbital computation b eing fully 3D. 



We use a recent determination of the absolute pro per motion of M 67 gi ven by Bellini et al. 



( 12010bl ). an d recent revisions of the Solar velocity ([Schonrich et al.ll2010f) and oth er Galactic 



properties ( iBenjamin et al.ll2005l ; IChurchwell et al. 



2009 



Brunthaler et al.ll2011f ). 



This paper is organized as follows. In Section [2] we provide data for the Sun and 
M 67 used to compute their Galactic orbits. The Galactic mass distribution used in our 
computations is presented in Section |3j The employed numerical scheme is given in Section 
HJ We present and discuss our results in Section |5j Our conclusions are given in Section |6j 



2. M 67-Sun parameters 



M 67 (NGC 2682) is a well studied, nearby (~900 pc) cluster. Its proximity has 



made it an ideal target for a wide range of studies from radial velocities (IMathieu et al. 



19861 ) and proper motions (lSanderslll977l; iGirard et al.lll989t IBellini et al.ll2010bl ). to late- 



type stellar evolution, binary fractions ([Montgomery et al.l Il993l ) and dynamical structure 
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( IDavenport fc Sandquistll2010l ). Its resemblance to the Sun in age and metallicity makes this 
cluster one of the key objects in our Galaxy. 



In a recent paper, iBellini et al.l ( 12 01 Obi ) made use of the state-of-the-art ground-based 



reduction techniques for ground-based wide-field imagers ( lAnderson et al.ll2006 



Yadav et al. 



20081 ; IBellini et al.l 120091 ) to determine for the first time the absolute proper motion of M 67 
using faint background galaxies as reference. In this work we decided to adopt their deter- 
mination, since it is a pure differential measurement and does not rely, as other works do, 
on complex r egistra tions to the ICRS system through a global network of objects. From 
Bellini et al.l ( )2010bl ). we list in Table [1] the absolute proper motion of M 67, and also its 
radial velocity, vy, and distance, r, with respect to the Sun, along with their corresponding 
uncertainties. 

In our computations, four parameters needed for the Sun are its velocity components 
with respect to the local standard of rest, (U, V, W)q, and its position perp endicular to the 
Galac tic plane, z . The Solar velocity (U, V, W) Q has recently been revised by lSchonrich et al. 
(120 lOl ). with a significant increase in the value of V Q , in the direction of Galactic rotation. 
This Solar velocity is listed in Table [TJ taking the U comp onent negative towards the Galac- 
tic center. The listed uncertainties are as considered by iBrunthaler et al.l (120111 ). who use 
the new Solar velocity to obtain a revised value for the ratio of 0o and Rq, the circular 
rotation speed and Galactocentric di stance of the loca l stan dard of rest. This ratio and a 
weighted average of Rq, also given by IBrunthaler et al.l (|2011l). l ead to a corresponding value 
of ©o- These values of Rq and Go given by IBrunthaler et al.l (120111 ) are also employed in 



our computations, a nd are considered in the fo llowing Section. They are within the wide 
intervals obtained by iMcMillan &: Binneyl (120101 ) from an analysis of motions of a sample of 
masers in star-forming regions. 



The z Q position is discussed by iReedl (120061 ). giving several recent estimates that have 
been obtained by different authors. From his Table 1 we have taken the representative 
z -interval as 20±5 pc. With a given initial (i.e. present) z position of the Sun in the 
considered Galactic potential, we add z Q to the initial z position of M 67, obtained with 
the usual convention that the Sun is on the Galactic plane. The z Q interval is also listed in 
Table [U 



3. The Galactic model 



We have employed a detailed Milky- Way semi-analytic model to compute the orbits 
of the Sun and M 67. This model has axisymmetric and non-axisymmetric components, 
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the latter due to spiral arms and a Gal actic bar. To build the Galactic model we start 



with the axisymmetric Galactic model of lAUen fc Santillanl ( Il99ll ). which consists of three 
components: disk, spherical bulge, and spherical halo. All the mass in the spherical bulge 
is now employed to build the Galactic bar, and a small fraction, discussed below, of the 
total mass of the disk is employed to build the spiral arms. Thus, the only axisymmetric 
components in the final model are the diminished disk and the spherical halo. 

The inclusion of the spiral arms and the Galactic bar does not modify the total mass 
of the original axisymmetric model. Thus, the mean circular rotation speed on the Galactic 
plane is maintained; in particular 0o ~ 220 km/s at a Galactocentr i c dist ance Rq = 8.5 
kpc, which are the values of these parameters in the lAllen fc Santillanl ( Il99ll ) model. These 
values of -Rn and 8 n ar e now modified in the non-axisymmetric model to those given by 



Brunthaler et al.l (120111): Ra = 8.3± 0.23 kpc, 9 = 239±7 km/s. This is done by first 
scaling the lAllen fc Santillanl (Il99ll ) model to a pair of values (i?o,@o) coming from the 



Monte Carlo sampling in our numerical scheme (discussed in Section H]), thus modifying the 
total masses of the original components, following with the transformation of the final bulge 
into the Galactic bar, and taking a fraction of the final disk to build the spiral arms. 

In the following two subsections we comment on some parameters of the Galactic spiral 
arms and the Galactic bar. Table [2] lists the values of these parameters. 



3.1. The Galactic Spiral Arms 



3.1.1. Geometry 



From the Spitzer /GLIMPSE database, Benjamin et al.l ( 120051 ) and IChurchwell et al. 



( 120091 ) have given new results on the structure of the Galactic s piral arms, as well as th e 
Galactic bars in the inner region of our Galaxy. As pointed out by IChurchwell et al.l ( 120091 ). 
although the geometry of the Galactic spiral structure is the most problematic, it appears 
that our Galaxy has two grand-design spiral arms: the Scutum-Centaurus and Perseus arms, 
associated with overdensities in the old stellar disk, and two additi onal secondary arms, the 
Sagitt arius and Norma arms, associated with gas overdensities. ISteiman-Cameron et al. 
( 20101) discuss the difference between these two types of arms, suggesting, as in lMartos et al. 
( 20041 ). that the gaseous spiral arms are the response of the gas to the main stellar arms. 



In our model we include only a spiral pattern that represents the two main stellar arms. 
Our three dimensional (3D) orbital compu tations require a 3D p otential for these arms. 
We adopted the spiral arm model given by iPichardo et al.l (120031 ). which consists of a 3D 
superposition of inhomogeneous oblate spheroids along a given spiral locus, adjustable to 
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better represent the available observations of the Galactic spiral arms. 



In Figured] we give details of the geometry of the spiral arms in our model. T his distance- 
scaled figure on the Galactic plane is a reproduction of figure 2 in iDrimmell (120001 ) . The spiral 
traces shown with black squares are the observed gaseous spiral arms. The two continuous 
lines, on which we illustrate t he superposition of spheroids (shown as circles), are the two 
stellar spiral arms obtained by IDrimmell (120001 ) in the K band, with a pitch angle of 15.5°. 
These stellar arms are not those stellar arms obtained in the Spitzer / G LI M P S E data; the 
latter correspond to the arms marked with 1 (Perseus) and 2 (Scutum-Centaurus) in Figure 
[TJ (thus, there are two gaseous arms that approximately coincide with the two main stellar 
arms). However, due to the uncertain geometry of the spiral structure (jChurchwell et al. 



20091 ). we have considered an interval of values for the pitch angle of our two spiral arms. 



In Figure [TJ the short-dashed lines mark a lower limit of 12° in the pitch angle, and the 
l ong-d ashed lines an upper limit of 19°, as in the second solution obtained by iDrimmel 
(120001 ) (but with a different spiral locus). Thus, in our numerical simulations we take the 
pitch angle of the spiral arms as i = 15.5±3.5°; the two Spitzer / G LI M P S E stellar arms are 
approximately contained in this pitch-angle range. 



The radial-extent interval of the spiral arms is taken as [2.6,12] kpc, with i?o 
This interval is scaled with the value of Rq given by the Monte Carlo sampling. 



8.5 kpc. 



3.1.2. Strength 



The strength of the spiral arms is related with their total mass, which is a small fraction 
of the disk's mass. The distribution of this mass along the arms i s taken as exponentia l, 
falling with the radial scale length of the exponential d isk modeled bvlBeniamin et al.l ( 120051 ): 



H* = 3.9±0.6 kpc , using R = 8.5 kpc (although the iMiyamoto &: Nagail ( 119751 ) disk in the 



Allen &: Santillanl (Il99lf ) model has not strictly an exponential density fall). 



To quantify the strength o f the arms we compute the function Q T ([Sanders fc Tubbs 



19801 ; ICombes fc Sanders! 1 198 11 1, which is the ratio of the maximum azimuthal force of the 
spiral arms at a given Galactocentric distance on the Galactic plane, to the radial axisym- 
metric force at that distance. The maximum value of Qt over the radial extent of the spiral 
arms, called Q s = (Qr)max; , is a measu r e of t he st rength of the spir al arms. This param- 
eter has been computed by iButa et al.l ( 12004 ) and iButa et al.l ( 120051 ) in a sample of spiral 
galaxies. 



To exemplify in our model the typical values expected for Q SJ we take the central values 
Rq = 8.3 kpc, O = 239 km/s of the (R ,Q ) distributions, and find Q s for an interval in 
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the ratio M arms /Mdi S k. Figure [2] shows the results for three values of the pitch angle i. 

In the ir sam ple of 147 spiral galaxies, which is dominated by Sbc and Sc galaxies, 



But a et al.l ( 120051 ) find that 75% of them have Q s < 0.20 . For Sbc galaxies, as our Galaxy, 
Q s is approximately less than 0.25 . Thus the intervals in M arms /Mdi S k and i considered in 
Figure [2] are appropriate for the Galactic spiral arms in our model. Specifically, we take 
Manna/Mdisk = 0.04±0.01; the Gaussian sampling considered in our computations (Section 
H]) will extend this range, and also the range in i, allowing for acceptable values of Q s in 



Figure[2J With this ratio M arms /Mdi S k and the la variations in R , © , the mass in the stellar 

©■ 



spiral arms lies in the interval 2.7 - 5.4 x 10 9 M e 



3.1.3. Angular velocity 



Gerhardl (120111 ) has given a recent review of different methods to determine the angular 
velocity of the Galactic spiral arms, Qs, as we U as those for the angular velocity of the 
Galactic bar. Some results are: Qs = 24 - 26 kms _1 kpc _1 finding the birt hplaces of open 



2003: Martos et al 



clust e rs (Dias & Lepine 20051 ) : Qs — 20, 40 kms *kpc 1 with gas flow models ( Bissantz et al. 



2004); Qs = 20 - 30 kms 1 kpc 1 with kinematics of OB and Cepheids 



stars ([Fernandez et al.ll200ll ): Qs ~ 18 kms _1 kpc _1 finding two families of closed orbits in 
the reference frame of the spiral arms , which could be associate d with the Pleiades/Hyades 
and Coma Berenices moving groups ( Quillen fc Minchevl 120051 ) : Qs = 18, 25 kms~ 1 k pc~ 1 
from kinematic fits to five moving groups in the Solar neighborhood (IChakrabartyll2007l ) . To 
represent all these results, in our numerical simulations we take Qs 
as listed in Table [2J 



24±6 kms _1 kp< 



3.2. The Galactic Bar 



3.2.1. Geometry and strength 



For the Galactic bar, we co nsider the pro l ate po tential given by lPichardo et al.l (120041 ). 
which approximates a model of iFreudenreichl (119981 ) of COBE/DIRBE observations of the 
Galactic center. In our computations the specific 3D shape of the bar, prolate or triaxial, is 
not very important, because the orbits of the Sun and of M 67 will lie far from the region 
of the bar. A prolate model is therefore convenient to facilitate the computations (also, the 
COBE/DIRBE Galactic bar approximates a prolate figure). In Table [2] we list the bar's 
semi-major axis, scale lengths, and axial ratio; the scale length of the prolate bar along 
its minor axis is taken as the mean of those along the two minor axes of the triaxial bar 
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models of iFreudenreichl (119981 ). All the lengths of the bar will change under the Monte Carlo 
sampling of Rq. The la variations in R , B give a total mass of the bar of 1.5 - 1.8 x 10 10 



MfTj, which lies in the range 1 - 2 x 10 10 of es t imated values (e.g..lKent 



Dwek et al.lll995l : lBlumlll995l : IStanek et al.lll997l : IWeiner fc Sellwoodlll999l ). 



1992: Zhao 1994 



3.2.2. Orientation and angular velocity 



The present orientation of the bar's major axis has been determined in several studies 



(e.g.JBinnev at all l991 



2005 



Weiner fc Sellwood 1999; Bissantz fc Gerhard 2002; Babusiaux fc Gilmore 



Minchev et al.l 120071 ). The average of the mean values of t he angle betwe en the bar's 
major axis and the Sun-Galactic center line is around ~ 25°. iGerhardl (120021 ) suggests </> 



= 20° as a good working value. We take this angle for the present-day orientation of the 
Galactic bar. 



There is also a long list of studies to estima te the bar's angular velocity, Qb, (IGerhard 



20 111 , and references therein). I Gerhard! (120111 ) concludes from his review that the most 
likely range in is 50 - 60 kms _1 kpc _1 . Thus, in our computations we take = 55±5 
kms _1 kpc _1 , as listed in Table [2j 



The numerical scheme 



The parameters listed in Tables [T] and [2J with their corresponding uncertainties, are 
employed in the numerical simulations. The listed uncertainties are considered as la vari- 
ations, and a Gaussian Monte Carlo sampling generates the parameters to compute the 
present-day positions and velocities of the Sun and M 67, along with the needed parameters 
in the Galactic potential. There are 14 variables to be sampled, 8 from Table (U and 6 
from Table [2J For a sampled Rq, the lengths associated with the bar were transformed, and 
the scale length of the spiral arms, H+, was taken as the corresponding scaled value of the 
mean H+ = 3.9 kpc wit h Rp = 8.5 kpc. W e used the Gaussian random-number-generator 



routine gasdev given by lPress et al.l (119921 ) . For each generated set of parameters the orbits 
of the Sun and M 67 were simultaneously computed backward in time, in the time-varying 
Galactic potential, up to the Solar age, 4.57 Gyr, allowi ng for a 0.2 Gyr u ncertainty. The 
high-precision Bulirsch-Stoer algorithm implemented by iPress et al.l ( 119921 ) was used in all 
the orbital computations. 

We ran N tot = 3.5 x 10 5 pairs of orbits Sun-M 67, searching for close encounters in the 
past; especially in the assigned time interval for the age of the Sun, (At)g un = — 4.57±0.2 
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Gyr. However, close encounters outside this interval were also considered relevant, since, 
as commented in Section 15.11 the Solar System could live within M 67 during 2 or 3 Gyr 
before the Sun could escape from this cluster, thus occurring a close encounter Sun-M 67 in 
a backward-time computation. 

Two major assumptions in our computations are that the structure of the spiral arms 
and the bar, as well as their angular velocities, remain unchanged during the 4.77 Gyr covered 
in each run. There are some studies which suggest that the Galactic bar may be a long- 



lived, few Gyr old structure (INg et al.lll996l : ISevensterill999l ; ICole & Weinberg! l2002f ). This 



(Fuji et al. 


2011; 


Fovle et al. 


2011; 


Sellwood 


2011) 



the spiral arms can produce a significant effect on the orbits of the Sun and M 67, and our 
results suggest that the initial hypothesis of a common origin of the Sun and M 67 might be 
more sustainable taking the strength of the spiral arms as time dependent; thus modeling 
these arms as transient features. 



5. Results and discussion 

5.1. The Sun as a sibling of M 67's stars ? 

If the Sun is a sibling of M 67's stars, a convenient criterion for close encounters Sun- 
M 67 along their orbits through the Galaxy, is that the minimum approach distance, d min , is 
less than the tidal radius of the cluster; i.e., at the time of the en counter the Sun is w ithin 



the region of bounded stellar motions of M 67's stars. Recently, iHurley et al.l (120051 ) have 
computed N-body models for M 67. Their best model 2, which reproduces many properties 
of the present state of M 67, has an initial tidal radius of 31.8 pc. During the cluster's 
evolution this tidal radius decreases, due to the interaction with the Galactic potential. 
Thus, to consider the early phases in the cluster's evolution, we analyze close encounters for 
which d m[n < 20 pc and d min < 30 pc. 

To illustrate the form of the resulting distribution of close encounters in our compu- 
tations, up to d min = 100 pc, Figure [3] shows the logarithm of the distribution of close 
encounters Sun-M 67, f(rf m in) = (number of encounters in the <i m i n interval) /N to t, N tot = 3.5 
x 10 5 . This figure includes all the encounters in a given d m i n interval, independently of the 
time they occurred. The continuous line is for 10-pc bins in d m \ n) the short-dashed line for 
20-pc bins, and the long-dashed line for 30-pc bins. As a start, there are many encounters 
in the intervals of interest, d min < 20 pc, d min < 30 pc. 

Taking from Figure |3] only the close encounters with d m[n < 20 pc, and d m i n < 30 pc, 
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we show in Figure H] their distribution in relative velocity between the Sun and M 67 at the 
time of the encounter, f(v Te \) = (number of encounters in the i> re i interval, with (i min < 20 
pc or d min < 30 pc)/N tot . The peaks in the distributions are around 50 - 60 km/s, and no 
encounters were obtained with u re i < 20 km/s. 

The same analysis can be made especially for those close encounters which occurred 
within the assumed time interval for the Sun's age, (At) Sun = — 4.57±0.2 Gyr. Figure [5] 
shows their distributions in relative velocity, with d min < 20 pc, and d min < 30 pc. The 
majority of these encounters have relative velocities larger than 30 km/s, and as already 
shown in Figure HI there are no encounters with v re \ < 20 km/s. 



With the assumption that the Sun was born in a star cluster, lAdamsl (120101 ) estimates 
that the planetary system, formed in a time scale about 10 Myr, can survive within the 
cluster for ~ 0.25 Gyr if the cluster's star density is n* ~ 10 3 pc ~ 3 . Th is surviving time 



increases if n+ decreases. The M 67 N-body model of iHurley et al.l (120051 ) has an initial n+ 
= 50 pc~ 3 within the half-mass radius, i.e. M 67 was not an initial dense system. With this 
n±, the surviving time for the Solar System within M 67 could be 2 or 3 Gyr, followed by the 
escape of the Sun. Thus, considering close encounters outside the Sun's-age time interval is 
relevant in this respect. However, the important constraint of the survival of the planetary 
system imposes a corresponding constraint on the relative velocity between the Sun and 
M 67 at the time of the Sun's ejection (encounter, in our backward-time computation). 



5.1.1. The Sun ejected from M 67 by a three-body encounter ? 



The limit v re \ > 20 km/s obtained in our computations results a high value compared 
to the smooth escape-velocity ejecti on from a star clus ter, with a typical internal velocity 
dispersion of ~ 1 km/s. According to lGirard et al.l ( 119891 ). the present-day velocity dispersion 



in M 67 is 0.81±0.10 km/s. However, in these cluster environments, three-body encounters 
are typical processes capable to produce high-velocity ejections; p articularly a binary-s i ngle 



1983a 


b, 


1984. 


1985. 


1993; 


Heggie et al. 


1996) 



another star in the cluster, a bound triple system is formed, and after some time one o f the 



stars is ejected (this process is known as resonance scattering; iHut fc Bahcall (1983)); or 
the captured star is ejected in the interaction (this is called a flyby; IHut fc Bahcalll ( 119831 ) ) . 
For high velocity ejections, hard binaries play the important role. In this case the ejected 
star leaves the three-body system with a velocity of the ord er of the intern al binary velocities 



( IHutl 119931 ). Thus, this ejection velocity is approximately (lAarsethl 120061 ) 
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'-'ejec — 



/G(mi 



+ m 2 , 



2a 



1/2 



(1) 



with mi, W2 the masses of the stars in the binary system, and a its semi- major axis. Thus, 
for a hard binary with equal mass components, and individual masses of order ~ 1 M , the 
ejection velocity is of order (GM Q / 'a) 1 / 2 . Taking this ejection velocity t> e j ec > 20 km/s, as 
obtained in our numer i cal sim ulations, it results approximately a < 2 AU (the M 67 N-body 
model of 
system) . 



Hurlev et al.l (|2005h h as an upper cut-off at a = 50 AU for the initial binaries in the 
Hut fc InagakiTl985 ) show that in an encounter of a hard binary with a single star, 



the typical minimum close distance between any two of the three stars can be less than the 
original semi-major axis a of the binary. In a resonance scattering this minimum distance 
can be of the order 0.01a. Thus, if the Sun were to be captured by a hard binary with rj a < 
2 AU, as obtained from our results, this small-distance interaction would destroy an initial 
circumstellar disk around the Sun, and no plane t s wou ld be formed; or the planets would 
be dispersed if they are already formed. lAdamsl (120101 ) comments that since the radius of 
the planetary system is ~ 30 AU (the position of Neptune), and a close encounter truncates 
an existing initial disk at a distance ~ 1/3 of its impact parameter, the required minimum 
distance in a close encounter suffered by the Solar System is ~ 90 AU; much larger than 
the < 2 AU in the triple-system interaction. Even with an order of magnitude increase in 
the total mass of the binary star, the semi-major axis, a < 20 AU, is still below the 90 AU 
limit. The situation would be worse if the ejected Sun from the triple system comes from a 
resonance scattering, where the strong interaction between th e three bodies may last orders 
of magnitud longer than the initial binary period (jHutlll993l ). and there are frequent close 
approaches. 

Thus, with the obtained results, we conclude that the Sun was not ejected from M 67 by 
a three-body encounter. Although we do find close encounters Sun-M 67 with an appropriate 
d m i n , the relative velocity of the encounter is too high. 



5.1.2. The Sun ejected from M 67 by an encounter with a giant molecular cloud ? 



A high star-ejection veloci ty from a star clus ter can also be produced b y an encounter 



with a giant molecular cloud (iGieles et al.l 120061 . and references therein) . ISolomon et al. 



(119791 ) find that the mass of the se clouds, M n , is in the range 10 5 - 3 x 10 6 M Q , and 
assuming their virial equilibrium ISolomon et al.l (119871 ) give a mass-radius relation M n = 
540i?^ M Q , with the clouds's radius R n in parsecs. 

Analytic estimations of the change in velocity of a star in the star cluster, at a distance 
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R* from its center, due to the encounter with a giant molecular cloud, can be obtained 
in the cases of distant and head-on encounters. For a distant encounter, in the impulse 
approximation, t he change in velocity in a direction perpendicular to the relative velocity 
cloud-cluster, is ( ISpitzerlll987l ) 



Av 



2GM n R+ 
b 2 V 



(2) 



with b the impact parameter, and V the magnitude of the relative velocity cloud-cluster. 
With Plummer models for bot h cloud and clust er, with median radii a\, 0,2, this relation 
applies with b > 5max(ai, 02) (iGieles et al.ll2006l ). 



For a head-on encounter, 6 = 0, and using Plummer models, the change in vel ocity in a 
direction perpendicular to the relative velocity cloud-cluster, is (IGieles et all 120061 ) 



Av 



2GM n R± 
V{Rl + al) 



(3) 



with V as in Eq. (j2|), and a n = R n /2. 



Taking the upper limit M n w 3 x 10 6 M for the mass of the giant mole cular cloud, the 
corresponding radius obtained with M n = 540i?^ M & ([Solomon et al.lll9871) is R„ ~ 75 pc . 
M 67 has an initial tidal radius of 31.8 pc, according to the model of iHurley et all ( 120051 ). 
Thus, Eq. (T5]) applies with b > 5a n = 5R n /2 « 190 pc. In our numerical simulations, at all 
close encounters Sun-M 67, which occur near the Galactic plane, M 67 has a velocity ~ 20 - 
40 km/s in the direction perpendicular to the Galactic plane (i.e., its P^-velocity component), 
and this is the dominant component of its peculiar velocity. On t he other hand, the typical 
random velocity of clouds is ~ 7 km/s (IBinney fc Tremaind 120081 ) . Thus, for V in Eqs. (j2j) 
and ()3]) we take V > 15 km/s. 

With b = 190 pc, V — 15 km/s, and .R* = 30 pc (i.e., a star, the Sun, at a distance 
from the center of M 67 equal to M 67's tidal radius), Eq. fl2]) gives Av < 1.4 km/s, and Eq. 
([3]), Av < 22 km/s. For encounters cloud-cluster other than distant or head-on, Av will lie 
between these estimates. 

Thus, a head-on encounter between M 67 and a most massive giant molecular cloud, 
with the Sun being a M 67's star at a distance equal to M 67's tidal radius, and located 
appropriately in a direction perpendicular to the relative velocity between M 67 and the 
cloud, could impart to the Sun an ejection velocity of order 20 km/s. These velocities are 
obtained in our computations, as shown in Figures H] and but with a low probability: less 
than 10 -4 for encounters Sun-M 67 occurring at any time (Figure H]), and less than 10~ 5 
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for encounters Sun-M 67 in the Sun's-age time interval (Figure [5]). The net probability is 
even lower, since all the above conditions on the type of encounter M 67-cloud, mass of the 
cloud, position of the Sun in M 67, and the Sun ejected appropriately for its orbit to be 
of low-z amplitude, must be fulfilled simultaneously. Furthermore, the a priori probability 
for the Sun being born in a star cluster with appropriate condition s consistent with present 
properties of the Solar System is 0.0085 (lAdams fc Laughlinl 120011 ). Thus, we conclude that 
there is a very low probability, much less than 10~ 7 , that the Sun was ejected from M 67 by 
an encounter with a giant molecular cloud. 



5.2. The Sun and M 67 born in the same molecular cloud ? 



As long as a molecular cloud has been assembled and no supernova e have been produced 
in its interior, the cloud may be approximately chemically homogeneous (iBland-Hawthorn et al 
20101 ). Thus, another possibility that can be analyzed with our computations is that the Sun 
and M 67 were born in the same molecular cloud. To test this possibility, we consider close- 
encounter distances Sun-M 67 up to 100 pc, say, and take only such encounters occurring 
within the time interval for the Sun's age, or subintervals centered at its mean age. Figure [6] 
shows their distributions in relative velocity between the Sun and M 67, for three time inter- 
vals centered at —4.57 Gyr: — 4.57±0.2 Gyr, continuous line; — 4.57±0.1 Gyr, short-dashed 
line; — 4.57±0.05 Gyr, long-dashed line. Only in the first interval there is a close encounter 
with v vc \ = 17.8 km/s; in all three intervals the rest of encounters have t> re i > 20 km/s. This 
is a high velocity compared with typical random velocities le ss than 10 km/s within clouds 
with sizes less than 100 pc ( lLarsonlll98ll ; ISolomon et al.lll987l ). Also, at all these encounters 
M 67 has a velocity ~ 20-40 km/s in the direction perpendicular to the Galactic plane, which 
also exceeds the typical 7 km/s random velocity between clouds. Thus, we also exclude the 
possibility that the Sun and M 67 were born in the same molecular cloud. 



5.3. The effect of the spiral arms 

The conclusions in Sections 15.11 and 15.21 depend strongly on the high relative velocity 
between the Sun and M 67 at their close encounters. The relative velocity, v re \, between 
the Sun and M 67 in a close encounter depends on their respective Galactic orbits. The 
high values of t> re i obtained in the computations reflect mainly the different z amplitudes 
above or below the Galactic plane reached by their orbits. At the present time, M 67 is 
above the Galactic plane, at z w 400 pc; the Sun is at z « 20 pc. If we compute their 
orbits in an axisymmetric potential (both will be box type), the Sun would stay within \z\ 
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~ 50 pc, and M 67 within |z| « 400 pc. Any close encounter obtained in this case will have 
a high v ie \, since M 67 will have a high z-velocity when it crosses the low-z region where 
the Sun moves. A non-axisymmetric Galactic potential is more promising to investigate the 
hypothesis of a common origin Sun-M 67, because the Galactic orbits are perturbed by the 
non-axisymmetric co mponents, and v rP i can reach lower values. For our problem, the spiral 
arms play a key role. lAUen et al.l (120081 ) have investigated the effect of the spiral arms, of 
the type we are using in the present study, on the orbits of some globular clusters. They 
used M arms /M disk = 0.03, and found a slight difference in z amplitudes compared with the 
orbits computed in an axisymmetric potential; the main difference appeared in the radial 
motions, parallel to the Galactic plane. However, in our present study we are extending the 
range of possible values of M arms /Mdi S k; values as high as 0.06 - 0.07 are being sampled by 
the Monte Carlo scheme, consistent with a Q s for a galaxy with the same Hubble type as 
the Milky Way (Figure EJ. 

Figure [7J shows the meridional orbits of the Sun (low-z amplitude) and M 67 (high-z 
amplitude) for one of our runs. This case has M arms /Mdi S k = 0.0123; i.e., a low strength of 
the spiral arms. The black dot shows the position of the close encounter. Note that there 
is a slight perturbation on both orbits. In this case the relative velocity at the encounter is 
u re i = 49.6 km/s; a high value. 

With the same (R,z) region as in Figure [TJ we show in Figure [8] another case, but now 
•Manng/Mdisk = 0.0636 . This is a case with high strength of the spiral arms. The present-day 
positions of the Sun and M 67 are marked with ®, and the black dot shows the corresponding 
position of the close encounter. Here, the orbit of the Sun shows a strong perturbation in 
the radial direction, and M 67's orbit has also strong perturbations in both radial and z 
directions. Note especially that the z amplitude, |z| max , in the orbit of M 67 around the time 
of the close encounter, is nearly half of that reached in Figure [7J In this case the relative 
velocity at the encounter is v Te \ = 24.3 km/s, one of the lowest values of v re \ obtained in our 
computations. 

If it were possible to find situations in which the z amplitude, |z| max , of M 67's orbit 
around the time of the close encounter, were similar to the z amplitude of the Sun, then 
v Te [ would decrease, and more favorable conditions might be found to support a common 
origin of the Sun and M 67. In this respect, an important question is: how has M 67 
acquired its present, high- z position? One possibl e answer is that it keeps memory of 
its birth at such high z. Vande Putte et al. ( 201oh a nalyzed a sample of open clusters 



from the open cluster catalogue by 
ters in the sample have l o w \z ' 



Dias et al.l (120021 ). The majority of the open clus- 
90% have |z| max < 350 pc. For clusters with higher 
Vande Putte et al.l (120101 ) mention the proposed formation mechanisms at high z 



max' 
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of iMartos et all ( 119991 ) and Ide la Fuente Marcos fc de la Fuente Marcos! (|2008h . and also 
suggest others. From the WEBDA Open Cluster Database (IPaunzen fc Mermilliodl 120081 ) 



de la Fuente Marcos fc de la Fuente Marcos! (120081 ) give tw o examples of young clusters that 



are likely born at high-z distance from the Galactic plane. iFriell ( 119951 ) also discusses high-z 
distance formation scenarios for the old open clusters, like M 67. Thus, M 67 could have 
been born at a |z| distance comparable with its present z « 400 pc. But in this case, if 
the Sun were to be ejected from M 67, a high t> rel would be needed to put it into its low-z 
amplitude orbit, and we obtain the conclusions in Sections 15. II and I5T21 

Another possible answer for the present, high-z, position of M 67 is that it was born in 
a low-z amplitude orbit, and has been dispersed in the z direction by the force of the spiral 
arms, and collisions with interstellar clouds. M 67 has been able to withs tand the tidal force 
of the clouds without being disrupted, as theory predicts (ISpitzerlll958l ). As we have said 
above, this situation might be more favorable to a common origin of the Sun and M 67; the 
Sun leaving M 67 with a lower value of v re \ when M 67 is in its low-z amplitude phase. 

The example presented in Figure El shows that there is an increase in the z amplitude of 
M 67's orbit from the time of the encounter to the present time; i.e., the orbit is dispersed to 
high z in this interval of time. Figure |9] shows the evolution of the coordinate z as a function 
of time in M 67's orbit for the case in Figure [HI but the time extended up to 5 Gyr in the 
past, to show more clearly the effect of the spiral arms. The black dot shows the coordinates 
(t,z) of the close encounter Sun-M 67, which occurred around 1 Gyr in the past. There is 
approximately a periodic increase- decrease in this z coordinate. The example in this Figure 
|8] and cases with similar behavior in the z coordinate, suggest that maybe considering non 
steady spiral arms, their increase- decrease periodicity in the z coordinate could be changed, 
obtaining situations in which this coordinate further decreases. This tr ansient-spiral-arm s 
mod eling wou l d be m ore in line with the conclusions of iFuji et al.l ( l201ll ). iFoyle et al.l ( 120111 ). 
and ISellwoodl (1201 lh . This will be done in a future study. 



6. Conclusions 

Considering the resemblance in age, metallicity, and distance from the Galactic center, 
of the old open cluster M 67 and the Sun, we explored in this paper the possibility that the 
Sun was once a member of M 67. We employed for this purpose the most recent proper- 
motion determination of M 67 in orbital computations in a 3D non-axisymmetric Milky Way 
model. The non-axisymmetric Galactic model includes a bar and 3D spiral arms, rotating 
with their particular angular velocities. The employed parameters to specify the properties 
of the bar and spiral arms have been taken from several papers from the literature. We used 
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a Gaussian Monte Carlo sampling to generate different present-day orbital conditions for the 
Sun and M 67, and parameters of the Galactic model. We ran 3.5 x 10 5 simultaneously- 
computed pairs of orbits Sun-M 67, looking for close encounters in the past. The analysis of 
these close encounters shows that the corresponding relative velocity between the Sun and 
M 67 is larger than 20 km/s. This velocity is too high; a three-body encounter within M 67, 
with the Sun being one of the three bodies, and giving this ejection velocity to the Sun, would 
destroy an initial circumstellar disk around the Sun, or disperse its already formed planets. 
Thus, the Sun was not ejected by a three-body encounter in M 67. Also, by analyzing a 
possible encounter of M 67 with a giant molecular cloud, we find a very low probability, 
much less than 10~ 7 , that the Sun was ejected from M 67 by such an encounter. The high 
values of the relative velocity also exclude the possibility that the Sun and M 67 were born 
in the same molecular cloud. We have illustrated the effect of the spiral arms on the Galactic 
orbits of the Sun and M 67. Modeling the spiral arms as transient features might prove to 
be compatible with the Sun-M 67 common-origin hypothesis. 

We thank the anonymous referee for a very careful review and several excellent sug- 
gestions that greatly improved this work. B.P. and E.M. thank support by grant PAPIIT 
IN110711-2. A.B. acknowledges the support by MIUR under program PRIN2007 (prot. 
20075TPK9). 
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Table 1. M 67-Sun parameters 



Parameter Value References 



M 67 



(fl a COS 5)2000.0 


— 9.6±1.1 mas yr 1 


1 


(^5)2000.0 


— 3.7±0.8 mas yr -1 


1 




33.78±0.18 km s" 1 


1 


r 


815±81.5 pc 


1 


Sun 




(-11.1±1.2,12.24±2.1,7.25±0.6) km s" 1 


2 


Z Q 


20±5 pc 


3 



References. - 1) Bellini et all ( 12010bh . 2) ISchonrich et all ( 120101 ). 
3) iReedl (bood ). 
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Table 2. Sun's Galactocentric distance, local circular rotation speed, and parameters of 

the non-axisymmetric Galactic components 



Parameter 


Value 


References 


Ro 


8.3±0.23 kpc 


1 


e 


239±7 km/s 


1 


Spiral Arms 


pitch angle (i) 


15.5±3.5° 


2 


scale length (H±) 


3.9±0.6 kpc (.Rq = 8.5 kpc) 


3 


M aTms /Mdisk 


0.04±0.01 




mass 


2.7-5.4 x 10 9 M 




pattern speed (fts) 


24±6 kms^kpc -1 


4 


Prolate Bar 


semi-major axis 


3.13 kpc (Rq = 8.5 kpc) 


5 


scale lengths 


1.7, 0.54 kpc (R = 8.5 kpc) 


5 


axial ratio 


0.54/1.7 




mass 


1.5 - 1.8 x 10 10 M 




angle between major axis 






and the Sun-GC line 


20° 


6 


pattern speed (fig) 


55±5 kms _1 kpc _1 


4 



3) 
6) 



'. References. 1} 



References. - 1) iBrunthaler et al.l (120111) . 2) iDrimmell (2000). 
Benjamin et -all d2005h. 4) iGerhardl feoilh . 5) iFreudenreichl Jl998[ ). 
Gerhard! tom . 



J I I L 



J I I L 



J L 







x/R 



Fig. 1. — Reproduction of figure 2 in IDrimmell ( l200ol ). illustrating in our model the superpo- 
sition of spheroids (circles) along t he stellar spiral l ocus (continuous lines, with pitch angle 
15.5°) in the K band, obtained by IDrimmell ( 120001 ) . The short-dashed lines mark a pitch 
angle limit of 12°, and the long-dashed lines a 19° limit. 
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0.02 0.04 0.06 

M arms/ M disk 

Fig. 2. — The parameter Q s , which gives a measure of the strength of the spiral arms, as a 
function of the ratio M arms /Mdi S k for three values of the pitch angle i. In this figure we take 
Ro = 8.3 kpc and 6 = 239 km/s. 



-29 - 



-1.5 - 



-2 - 



c 

• i— i 



-2.5 - 



O 



-3 - 



-3.5 - 








20 



40 60 
d min (P C ) 



80 



100 



Fig. 3. — Logarithm of the close-encounters distribution f(d m i n ), for d m[n bins of 10 pc 
(continuous line), 20 pc (short-dashed line), and 30 pc (long-dashed line). The encounters 
occur at any time, up to the Sun's age. 
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Fig. 4. — The distributions in relative velocity between the Sun and M 67 for the close 
encounters in Figure [3] with d min < 20 pc and d m in < 30 pc. 
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Fig. 5. — The distributions in relative velocity between the Sun and M 67 for close encounters 
which occurred in the Sun's-age time interval (Ai)s un = — 4.57±0.2 Gyr. Encounters with 
^min < 20 pc and d min < 30 pc are considered. 
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Fig. 6. — Distributions in relative velocity between the Sun and M 67 for close encoun- 
ters with d m i n < 100 pc, occurring in the time intervals: — 4.57±0.2 Gyr, continuous line; 
— 4.57±0.1 Gyr, short-dashed line; — 4.57±0.05 Gyr, long-dashed line. 
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Fig. 7. — Meridional orbits of the Sun (low-z amplitude) and M 67 (high-z amplitude) for a 
run with a low strength of the spiral arms, M arms /Mdi S k = 0.0123. The black dot shows the 
position of the close encounter. 
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Fig. 8. — As in Figure [7J with the same (R,z) region, here for a run with a high strength of 
the spiral arms, M arms /Mdi S k = 0.0636 . The points <g> mark the present-day positions of the 
Sun and M 67, and the black dot shows the position of the close encounter. 
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Fig. 9. — Evolution of the coordinate z in the orbit of M 67, for the case shown in Figure EJ 
extending the time up to —5 Gyr. The black dot shows (t,z) at the close encounter Sun-M 67. 



